Correlation structure in nondipole photoionization 



co 
o 
o 

(N 
C 

a 



Si 



O 
43 

Oh: 



> 

o 
o 

o 

CO 

O 

O 

43 : 

Ph. 



M. 



Ya. Amusia, 1 * 2 A. S. Baltenkov, 3 L. V. Chernysheva, 2 Z. Felfli, 4 S. T. Manson, 5 and A. Z. Msezane 4 

1 Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 
2 A. F. Ioffe Physical- Technical Institute, St. Petersburg 194021, Russia 
3 U. A. Arifov Institute of Electronics, Tashkent 700187, Uzbekistan 
■* Center for Theoretical Studies of Physical Systems, Clark Atlanta University, Atlanta, GA 30314 
'Department of Physics and Astronomy, Georgia State University, Atlanta, GA 30303 

(Dated: February 2, 2008) 

The nondipole parameters that characterize the angular disribution of the photoelectrons from 
the 3d subshell of Cs are found to be altered qualitatively by the inclusion of correlation in the 
form of interchannel coupling between the 3d 3 / 2 and 3d 5 /2 photoionization channels. A prominent 
characteristic maximum is predicted only in the parameters for 3d 5 /2 photoionization, while the 
effect for 3d 3/ <2 is rather weak. The results are obtained within the framework of the Generalized 
Random Phase Approximation with Exchange (GRPAE), which in addition to the RPAE effects 
takes into account the rearrangement of all atomic electrons due to the creation of a 3d vacancy. 
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Nondipole effects in the photoionization of atoms and 
molecules were thought to be of any importance only 
at multi-keV photon energies[l, 2] despite indications to 
the contrary [3, 4, 5]. However, in an upsurge of work 
on the subject, it was found both theoretically and ex- 
perimentally that nondipole effects are often of signifi- 
cance at photon energies of hundreds and even tens of 
eV[6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]. It has been 
determined that photoelectron angular distributions are 
particularly affected by these nondipole photoionization 
channels. In lowest order, the effect of the non-dipole 
channels is to add two new terms to the well-known 
dipole expression for the photoelectron angular distri- 
bution, along with their associated dynamical parame- 
ters, 7 and 5, the nondipole parameters[5, 7, 8, 17, 18]. 
These parameters arise from the interference of dipole 
and quadrupole (El and E2) channels in much the same 
manner that the f3 parameter, which characterizes the 
dipole photoelectron angular distribution, arises from 
dipole-dipole interference[19]. 

Studies of these nondipole parameters have found 7 
and S to be essentially structureless, except in regions 
of dipole or quadrupole resonances [11, 16] or Cooper 
minima[13]. It has also been averred that correlation 
effects are unimportant for nondipole effects[20]. In this 
communication, we show that assertion to be incorrect. 
In fact, we report on a new phenomenon, a new kind of 
structure in nondipole parameters that has nothing to 
do with resonances or Cooper minima. This structure, 
as shall be shown below, is due solely to many-body cor- 
relations; it has no analogue in a single-particle picture. 
Furthermore, the structure is induced by relativistic in- 
teractions; it does not arise in a non-relativistic calcula- 
tion. Thus, the study of this structure is a direct measure 
of relativistic correlation effects. 

To understand the basic origin of this structure, note 



that it was demonstrated recently [21] that the interac- 
tion between photoionization channels belonging to dif- 
ferent components of the spin-orbit doublets 3c? 3 / 2 and 
3c?5/2 in Xe, Cs and Ba affect dramatically the partial 
photoionization cross sections. Specifically, it was shown 
that due to this interchannel coupling interaction, the 
partial 3c? 5 / 2 cross section acquires an additional promi- 
nent maximum. This result gave the explanation for the 
recent experimental observation of this effect in Xe[22] 
and also predicted similar, even more dramatic effects 
for 3c? photoionization in Cs and Ba. 

In this communication we present the results of our 
studies, based on the approach developed earlier[21], of 
the nondipole angular distribution parameters for 3c? 
photoionization in Cs to illustrate a new phenomenon 
where impressive manifestations of intra-doublet in- 
terchannel interactions producing new structures in 
nondipole parameters are predicted, thereby also con- 
tradicting the notion that correlation does not affect 
nondipole parameters significantly. 

The theoretical results were obtained using a spe- 
cially modified random-phase-approximation with ex- 
change (RPAE) approach previously developed for half- 
filled subshells[23, 24]. This calculation showed clearly 
the physics of the phenomenon observed [22] and led to 
results [21] that are in both qaulitative and quantita- 
tive agreement with experimental 3c? Xe data. Explic- 
itly relativistic calculations [25] were subsequently per- 
formed using the relativistic-random-phase approxima- 
tion (RRPA)[26] which confirmed the qualitative and 
quantitative accuracy of the results in [21]. 

To summarize the theoretical approach, we consider 
the 3c?5/2 and 3^3/2 subshells each as half- filled atomic 
subshclls. This permits us to apply straightforwardly the 
RPAE methodology to include many-body correlations 
(including interchannel coupling) for half-filled subshells. 
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Exchange is neglected between these two sorts of elec- 
trons, the six that form the 3d$/ 2 (called "up"), and the 
four forming 3d 3 / 2 (called "down") electrons. However, 
in the real half-filled 3d subshell one would have five elec- 
trons. But these corrections, 6/5 and 4/5, respectively 
can be introduced easily into the calculational scheme. 

Then we concentrate on the investigation of the influ- 
ence of " up" and " down" electrons upon each other and 
demonstrate that the effect of the 3d 3 / 2 ("down") elec- 
trons upon the 3c?5 /2 (" up" ) manifests itself not only in 
the partial cross sections ct 5 / 2 (uj), cr 3 /2(w)[21] and some- 
what in the respective angular anisotropy parameters 
(35/2(^)1 ^3/2 C'*') [25], but also in the nondipole parame- 
ters of the photoelectron angular distributions, 7(w) and 
S(uj). In all these cases the effect of the 3^/2 photoion- 
ization channels upon the 3d 5 / 2 leads to the creation of 
an additional maximum, while the action of the 3^/2 
electrons upon 3d 3 / 2 proved to be generally negligible. 
Note that since the matrix elements for the photoioniza- 
tion process are strongly modified by correlation in the 
form of interchannel coupling, it is evident that charac- 
teristics of the process other than the integrated cross 
section; for example, the photoelectron spin polarization 
(see [27] and references therein) are also modified. 

The angular distribution of photoclcctrons from an nl 
subshell created by linearly polarized light, including the 
lowest order nondipole parameters[5, 7, 8, 17, 18] is given 
by 



quadrupolc, qi± 2 ,o, arc defined as 
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where a n i{oj) is the partial cross section, f3 n i{uj) is the 
dipole angular anisotropy parameter, 7„z(w) and 8 n i(oj) 
are nondipole parameters, P 2 (cos6) is a Legendre poly- 
nomial, 9 is the angle between the photoelectron and the 
photon polarization directions, and <j) IS the angle be- 
tween the photon momentum and the projection of the 
photoelectron momentum in the plane perpendicular to 
the photon polarization. The expressions for 7„;(w) and 
8 n i(u>) in terms of dipole and quadrupole matrix elements 
were first obtained in a slightly different form in [5] and 
later in [7, 8, 13, 28]. Most of the attention in the past 
was given to I = and 1 = 1, where the expressions for 
-fni and S n i are relatively simple. Here we are interested 
in the more complex case of I — 2, and the explicit ex- 
pressions are given elsewhere [8]. For present purposes, 
however, it is necessary to point out that the expressions 
involve ratios of quadrupolc to dipole matrix elements, 
along with cosines of phase shift differences. 

In the present calculation, the 3d — > ep, ef dipole 
amplitudes and 3d — > es,ed,eg quadrupole amplitudes 
are taken into account for both 3d^/ 2 and 3d 3 / 2 pho- 
toionization. The matrix elements for dipole, dj±i, and 
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cfi±i = d n iM ±1 = f <f>ni{r)r<t> e i±i(r)dr, 
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$±2,0 = qnl,el±2,0 = J <l>nl{r)r 2 (j> € i±2fi{r)dr. (2) 



Here (f> n i (r) and (j)^ (r) are the radial parts of the Hartree- 
Fock (HF) one-electron wave functions. In the uncor- 
rected (HF) approximation, the formulae for 7 n ;(w) 
and 6 n i(u>) can be used directly with the dipole and 
quadrupolc amplitudes. 

An effect of correlation, however, is to render these 
transition matrix elements complex; the correlated 
dipole amplitudes are termed Di±\ and the correlated 
quadrupolc amplitudes Q;±2,o- To obtain the corre- 
sponding expressions for 7„ ; and S„i, one has to perform 
the following substitutions [5, 17] 

I d l±1 | 2 ->| ReD l±1 | 2 + I ImD l±1 | 2 
I ft±2,o | 2 H ReQi± 2fi | 2 + I ImQi ±2 , | 2 
rf/±i9/±2,ocos(<5 ;±2i o - 8 l±1 ) -» [ReD l±1 ReQi ±2 .o 
+ImDi ±1 ImQi ±2 .o)} cos(<S ;±2 ,o - 5j±i) - 
-[ReDi ±1 ImQi± 2 fi - ImDi ±1 ReQi± 2fi )} 

sin((5 ;±2:0 - S i±1 ) (3) 

where the 8v are the single-particle (HF) phase shifts for 
the designated channel, and Re and Im stand for real 
and imaginary parts, respectively. 

The method that we use here, just as in [21], is 
the spin-polarized random-phase approximation with ex- 
change (SPRPAE). However, for the intermediate 3d sub- 
shell SPRPAE is not sufficient; the effects of core rear- 
rangement (relaxation) must be taken into account. This 
is done by going from RPAE to the generalized (GRPAE) 
or, in our case, from SPRPAE to SPGRPAE which takes 
into account that while a slow photoelectron leaves the 
atom, the field seen is modified due to the alteration (re- 
laxation) of the wave functions of all other atomic elec- 
trons as a result of the creation of the inner-subshell va- 
cancy. GRPAE is discussed at length in [29]. Its exten- 
sion for a system with two types of electrons "up" and 
"down", the transition from SPRPAE to SPGRPAE, is 
straightforward. 

The results of our calculations for 7 for photoionization 
of Cs 3d 5 / 2 and 3d 3 / 2 are shown in Fig. 1; both the corre- 
lated and the uncorrelated Hartree-Fock (HF) results are 
shown to emphasize the features brought about by cor- 
relation. As seen from the figure, the outstanding differ- 
ence between correlated and uncorrelated results appears 
in the value of 7 for 3c? 5 / 2 photoionization in the region 
just above a photon energy of 740 eV. The fact that this 
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significant structure in 7 is not seen in the uncorrelated, 
HF, calculation is proof that its existence is due to cor- 
relation. Specifically, we have found that it is due to 
interchannel coupling among the 3cZ 5 / 2 and 3d 3 / 2 pho- 
toionization channels; by interchannel coupling we mean 
simply configuration interaction in the continnum. Since 
the " real" wave function for the final continuum state in 
a photoionization process is multichannel, there is mixed 
in with the final state wave functions for 3d 5 / 2 photoion- 
ization a small amount of 3<i 3 /2 channels. Then, due to 
the 3d — ► ef shape resonance in the region just above the 
3d 3 / 2 threshold, the 3d 3 / 2 dipole cross section is much 
larger than the 3d§/ 2 . Thus mixing a small amount of 
the wave function of the 3d 3 / 2 — > ef channel with the 
wave function of the much smaller 3d 5 / 2 channel alters 
the 3<i 5 / 2 dipole matrix elements significantly; the 3^/2 
cross section maximizes with the 3d 5 / 2 — ► ef just above 
its threshold, and drops off considerably at the energy of 
the 3d 3 / 2 — ► ef shape resonance. 

But the 7 parameter, as discussed above, is essentially 
a ratio of quadrupole to dipole matrix elements. There is 
also interchannel coupling among the quadrupole chan- 
nels, 3d 5 / 2 — > es,d,f and 3d 3 / 2 — > es,d,f. But there 
are no significant resonances in any of these channels in 
the energy range of interest, and all of the channels are 
of more or less the same size. Thus, interchannel cou- 
pling does not introduce any important changes in the 
quadrupole photoionization matrix elements like it does 
in the dipole case. Since only one participant in the ratio 
that yields 7 changes significantly because of correlation, 
then it is evident that the ratio itself must be altered 
as well. This is the reason for the significant structure 
seen in Fig. 1 for 7 of 3rf 5 / 2 in the vicinity of the 3d 3 / 2 
threshold. 

The only other real effect of correlation on the 7 pa- 
rameter is for 3d 3 /2 right at threshold, where the HF 
result is noticeably lower than the correlated value. This 
again is due to interchannel coupling, but in this case the 
effect of the 3g?5/2 photoionization channels on the 3c? 3 / 2 
causes the modification. In particular, the 3d 3 / 2 — > ef 
amplitude, which is extremely small at the 3d 3 / 2 thresh- 
old, is strongly modified by its interchannel interaction 
with the much larger 3^/2 amplitudes. 

The situation for the 5 parameter is shown in Fig. 2 
and we find here a story similar to that of the 7 pa- 



rameter. A significant correlation structure in the 3d 5 / 2 
channel, seen just above the 3d 3 / 2 threshold is evident, 
along with changes in 5 in the 3d 3 / 2 channel just at its 
threshold. These effects of correlation occur for exactly 
the same reason as discussed above for 7. It is apparent 
from comparison of Figs. 1 and 2 that the changes in 
5 due to correlation are qualitatively larger than for 7. 
This is probably due to the fact that the values of S arc 
roughly an order of magnitude smaller than 7. 

Owing to the geometry of a number of experimental 
setups to measure nondipole photoelectron angular dis- 
tributions, 7 and S are often not measured individually [9, 
10]; instead, ( = 7 + 3d is investigated. To connect 
with such experimental investigations, £ for the 3c?5/2 
and 3g? 3 /2 channels of Cs are shown in Fig. 3. Since £ is 
simply a linear combination of 7 and 5, the explanation 
for the correlation structure here simply follows from the 
previous discussion for 7 and S individually. 

In conclusion then, it has been shown that, due to spin- 
orbit induced interchannel coupling between the dipole 
photoionization channels arising from the two members 
of the 3d spin-orbit doublet in Cs, a new structure results 
in the nondipole photoelectron angular distribution pa- 
rameters 7 and S (and, of course, £). Similar results are 
found for Xc 3d and Ba 3d: these will be presented in a fu- 
ture publication. This is the first case of structure in the 
energy dependence of nondipole parameters discovered 
that is not related to resonances (dipole or quadrupole) 
or Cooper minima. And, it was shown explicitly that 
correlation in the form of interchannel coupling can in- 
deed be important for the nondipole photoelectron angu- 
lar distribution parameters, in contradistinction to the 
conventional wisdom. Experimental scrutiny into this 
prediction is highly desirable. 
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Figure Captions 

FIG. 1 Calculated values of the nondipole photoelectron 
angular distribution parameter 7 for Cs 3d 5 / 2 and 
3c? 3 /2 subshclls in correlated (SRPAE) and uncor- 
rected (HF) approximations. 

FIG. 2 Calculated values of the nondipole photoelectron 
angular distribution parameter 5 for Cs 3d 5 / 2 and 
3d 3 /2 subshclls in correlated (SRPAE) and uncor- 
rected (HF) approximations. 

FIG. 3 Calculated values of the nondipole photoelectron 
angular distribution parameter ( for Cs 3d 5 / 2 and 
3d 3 / 2 subshclls in correlated (SRPAE) and uncor- 
rected (HF) approximations. 
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